Dynamics of Cryogenic Jets: 
Non-Rayleigh Breakup and Onset of Non-Axisymmetric Motions 
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We report development of generators for periodic, satellite- free fluxes of mono-disperse drops with 
diameters down to 10 /mi from cryogenic liquids like H2, N2, Ar and Xe (and, as reference fluid, 
water). While the breakup of water jets can well be described by Rayleigh's linear theory, we find 
jet regimes for H2 and N2 which reveal deviations from this behavior. Thus, Rayleigh's theory 
is inappropriate for thin jets that exchange energy and/or mass with the surrounding medium. 
Moreover, at high evaporation rates, axial symmetry of the dynamics is lost. When the drops pass 
into vacuum, frozen pellets form due to surface evaporation. The narrow width of the pellet flux 
paves the way towards various industrial and scientific applications. 

PACS numbers: 07.20.Mc, 07.90,+c, 29.25.Pj, 47.55.Dz, 52.50.Dg 



Drop formation from incompressible liquids is ubiq- 
uitous in daily life and technology, and such processes 
have been at the focus of scientific investigations for al- 
most two centuries [l], Q, H, Q, H, @]. Even the simple 
case of jets emanating from a vibrating nozzle leads to 
a broad spectrum of drop sizes, and the existence of 
satellite drops [7]. However, for technologies like ink-jet 
printing 0] , for laser-plasma UV sources , for compact 
laser-based particle accelerators [Io|, accelerator exper- 
iments or for space operations [l2j], the drop sizes 
and rates must be highly homogeneous. 

The disintegration of a slow cylindrical jet is caused 
by the growth of perturbations initiated at the ejection 
nozzle. An axially symmetric perturbation leads to jet 
breakup when its amplitude becomes equal to the jet ra- 
dius. This domain is commonly denoted as the Rayleigh 
regime. In 1878 Lord Rayleigh described fluids in the 
limit of zero viscosity [2[, while his 1892 paper [3] con- 
cerns the opposite limit where inertial effects are neg- 
ligible compared to viscous ones. Rayleigh's work was 
extended by Chandrasekhar [4[ to arbitrary viscosities. 
It has been shown [7] that the measurement of breakup 
lengths, e.g. for the forced breakup of water jets with 
precise control of the nozzle-perturbation frequency, al- 
lows for precise tests of linear theories, although the final 
stage of capillary pinching is non-linear [l3j |. 

Experimentally, the selection of a well defined initial 
surface-perturbation amplitude and a certain frequency 
is commonly achieved by using piezoelectric transducers. 
Most of the existing data have been obtained for water 
(or water-based ink) and fluids of higher viscosities, with 
jets emanating into air at normal pressure. 

In this letter we report on first quantitative measure- 
ments of the breakup of cryogenic and water jets, in- 
jected into the same gaseous substance close to triple- 
point (TP) conditions (63 K and 125 mbar for N2, 14 K 



and 70 mbar for H2). In doing so we can widen the exper- 
imentally accessible parameter space, while keeping the 
Reynolds number Re=(p • Uj e t • Ro)/p (which represents 
the relative importance of inertial and viscous forces) and 
the Ohnesorge number Oh=p/ (a-p-Ro) 1 / 2 (viscous forces 
and surface tension) in the same range as for the existing 
water data, p (/i, a) denote the fluid density (viscosity, 
surface tension), and Vj et (Ro) the jet velocity (initial jet 
radius). At TP pressures and jet velocities of a few m/s, 
aerodynamic interactions with the ambient medium are 
small and, thus, we are close to the ideal case of an in- 
viscid jet breaking up in vacuum and our data allow for 
clean tests of Rayleigh's theory. 

According to the linear theory, a small surface pertur- 
bation 5 (measured in units of Rq) of a jet moving into 
z-direction, imposed at the nozzle exit (z=0) with fre- 
quency /, grows exponentially towards the jet axis 
Then, the minimum (over time) necking of the jet R{z) 
at a certain distance z from the nozzle is 
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The dimensionless growth rate 7/70 depends on the 
wavelength A of the perturbation and can be expressed 
in terms of a reduced wave number X=2ir • Rq/X = 
2tt • iV //vjet for X < 1 and Re • Oh > 1: 
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with 70 =0 h-(a/pRo) . Figure [T] shows a plot of Eq.© for 
water, N2 and H2 jets; linear theory predicts a practically 
identical X dependence of 7/70 for all three liquids. If Re- 
Oh approaches unity, the theory tends to underestimate 
7/70, see e.g. the discussion by Kalaaji et al. 

Since a jet breaks into drops for i?(z)=0, the jet length 
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FIG. 1: Dimensionless growth rate 7/70 as a function of the 
reduced wave number X in linear theory. For N2 and H2 two 
(practically identical) curves have been calculated with Eq. J2} 
using the Oh values from Table |U For water a typical Ohne- 
sorge number Oh=0.02 has been chosen, since 7/70 behaves 
similarly at Oh < 0.1. The arrows indicate the X- values for 
which we have observed satellite- free drop production. 



Lj e t results from Eq.([T]) as: 
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Relations ([I])-® have been experimentally verified in 
numerous studies with water or viscous fluids @, EEf. 
Satellite-free and mono-disperse water-drop production 
has also been reported: for not too large values of S 
(< 0.01) it onl y o ccurs at maximum values of 7/70, i-e. 
at X max =0.69 

Quantitative studies of forced jet break-up depend cru- 
cially on the suppression of unwanted nozzle vibrations 
which can, e.g., be caused by cold head units 0. Thus 
our drop generator [19[ is surrounded by a cryostat with 
baths of liquid N2 and He, see Fig. O Liquefaction of H2 
is achieved in three stages: first cooling with liquid N2, 
further cooling in a heat exchanger by evaporated He, 
and final cooling in the condenser by cold He gas. For 
production of, e.g., N2 and Ar jets, cooling with liquid 
N2 is sufficient while for water we use a dedicated genera- 
tor without cooling cryostat. The temperatures along the 
main gas and liquid flows are controlled with an accuracy 
of 0.1 K. From the condenser, the liquid passes through 
the vibrating nozzle and reaches the triple-point chamber 
(TPC) where the jet disintegrates into drops. The liquid 
temperature in the nozzle is maintained slightly below 
the boiling point and, according to simulation calcula- 
tions, the temperature decrease in the jet until break- 
up from heat exchange in the TPC is below 0.2 K (cf. 
Table [Tj). Therefore, the generation of semi-solid do- 
mains can be excluded. The temperatures and pressures 
in the TPC are kept over several hours within ±0.2 K 
(±0.2 mbar) at the nominal values. 

Two nozzle types have been utilized, glass nozzles in 
brass housings and such made from stainless steel. The 
former, with inner diameters of 12-30 jam at the nozzle 
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FIG. 2: (Color online) Sketch of the drop generator. The 
cryostat (total height ^0.8 m, not drawn to scale) comprises 
the cooling liquids (N2 and He) as well as the heat exchanger 
and the upper part of the condenser. These are cooled by 
a stream of evaporated He gas. The condenser houses a few 
cm 3 of the cryogenic liquid which is driven through a nozzle 
into the TPC by overpressure. Constant pressure in the TPC 
is maintained by an auxiliary gas feeding. 



tip, have the advantage of a smooth internal surface and 
allow one to look inside the channel during operation. 
The 16-30 /am steel nozzles offer high shape reproducibil- 
ity and smaller length-to-diameter ratios of the holes, 
allowing operation with lower jet-driving pressures. The 
results presented here have been obtained with glass noz- 
zles and pressures of ~0.4-0.9 bar for H2 and ~1.0-1.5 
bar for N2. Higher pressures lead to higher jet veloc- 



TABLE I: Parameters of the jets from Fig. [3] 7/70 is calcu- 
lated with Eq.©, the jet length L^ yleigh with Eq.®. 





N 2 


H 2 


Density p (kg/m a ) [20, 21] 


823.7 


64.04 


Static surface tension a (N/m) [20, 211 


0.0094 


0.0025 


Viscosity fi (mPa s) [20, 21] 


0.2113 


0.0183 


TPC pressure (mbar) 


300 


130 


TPC temperature (K) 


74 


17 


Jet temperature (K) 


77.2-77.0 


20.0-19.8 


Nozzle frequency / (kHz) 


26 


38 


Jet diameter 2Ro (/mi) 


17 


12 


Jet velocity Vj e t (m/s) 


2.6 


2.4 


Reduced wave number X 


0.53 


0.60 


Reynolds number Re 


84 


50 


Ohnesorge number Oh 


0.026 


0.019 


Re-Oh 


2.2 


1.0 


Growth rate 7/70 


0.31 


0.34 


Initial perturbation 5 


0.02 


0.02 


Expected jet length ^^ lei s h (^m) 


240 


70 


Measured jet length Lj e t (/am) 


310 


290 
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ities. The piezo-electric transducer allows us to excite 
sinusoidal nozzle vibrations in the range /=l-60kHz. 

The jet decay has been observed with two perpendic- 
ularly placed CCD cameras and strobe lamps that are 
triggered synchronously to the nozzle, with a 10 times 
lower frequency, and a flash duration of ~1.5 /is. The jet 
diameter 2Ro has been assumed to be equal to the nozzle 
diameter which is correct to better than ±10% {z|. This 
is confirmed by our measurement of jet diameters with 
the CCD cameras. Vj et has been determined by measur- 
ing the wave-propagation velocity along the jet surface. 
The initial amplitude 5 has been derived from reference 
measurements of water jet lengths using Eq.([3|), as well 
as from direct amplitude measurements with an Michel- 
son interferometer. See Ref. [15[ for further details of the 
jet diagnostics. 

Figure [3] shows the breakup of N2 and H2 jets; for 
certain choices of the jet and TPC parameters (see Ta- 
ble [Tj), satellite-free and mono-disperse drop production 
is achieved. Astonishingly, we find it at X=0.53 and 
X=0.60, respectively, i.e. for different wave numbers as 
compared to water. Since this is off the maximum from 
Fig.[lJ it suggests an X dependence of 7/70 different from 
Eq.([2|). In order to further test a possible deviation from 
Rayleigh-like behavior, we have checked the validity of 
Eqs. (pQ)-([3]) in more detail. 

The last two lines of Table U compare the measured 
and predicted jet-decay lengths; it is seen that the cryo- 
genic jets are much longer than expected. This cannot 
be explained by a variation of the parameters that en- 
ter Eq.(|3]). For example, an unphysical value of J=10 -8 
would be required to yield the large decay length of the 
H2 jet. The dynamic surface tension can differ from the 
static values for a quoted in Table H it had to be, how- 
ever, 20 (2) times smaller to explain the observed H2 (N2) 
jet lengths which seems unreasonable. 





« 



An explanation of the obvious violation of Rayleigh's 
theory might be disturbances like temperature, pressure 
and surface-tension variations as well as liquid-gas tran- 
sitions. It has, e.g., been predicted [22[ that evaporation 
and condensation decrease Lj et (which seems to contra- 
dict the large measured Lj et values from Table [Ij) and, 
for large evaporation rates, satellite production can be 
suppressed. One may presume that for our jets evapo- 
ration is strong since in the TPC jets and ambient gas 
coexist close to the vapor-pressure curve. In order to test 
the effect of evaporation we have enhanced it by reducing 
the jet velocities or the TPC pressure (keeping the other 
parameters constant). This leads to a phenomenon, to 
our knowledge not reported until now. The jets slowly 
move from their vertical flux direction to one off the noz- 
zle symmetry axis while preserving their smooth surfaces. 
The direction they choose seems to be random. However, 
once a jet has arrived at a certain stable configuration it 
may stay there for seconds; only at very small jet-bending 
radii are frequent abrupt changes of the jet directions ob- 
served. As an example we show in Fig. [H N2, H2 and 
water jets. Even back-bending is observed in some cases. 
Systematic studies reveal larger jet deflections from the 
vertical axis with decreasing jet diameter and velocity 
or TPC pressure. We note that asymmetric heating has 
been employed by another group to deflect liquid micro- 
jets |23[, however, we report spontaneous deflection of 



symmetrically produced jets. 

When the drops leave the TPC through a 1 st sluice into 
a subsequent chamber (p=o(10~ 2 ) mbar), they freeze to 
pellets due to strong surface evaporation, and are accel- 
erated by the gas flow from the TPC. Then the pellets 
pass a 2 nd sluice and a 2 nd chamber (o(10 -4 ) mbar), and 
finally reach a dummy scattering chamber through a thin 
tube (0^2 cm), which is located 1.2 m downstream of 
the TPC. In order to minimize turbulences of the gas 
flow in the 1 st sluice (and, thus, distortions of the pro- 
duced pellets) it has a circular cross section with a radius 
that decreases exponentially to 600 fim in flight direction. 
For pellet observation two CCD cameras have been po- 
sitioned at the outlet of the 1 st sluice or at the dummy 
chamber. Stable H2 and N2 pellet production with diam- 
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FIG. 3: First observation of satellite- free and monodisperse 
disintegration of N2 (left) and H2 (right) jets. In the upper 
edge of the photos the tip of the vibrating nozzle can be seen. 



2R Q = 17 |jm 
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1/^ = 4.0/3.9/3.0 m/s. 
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FIG. 4: First observation of non-axisymmetric N2 (a), H2 (b) 
and water (c-e) jets. 
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eters 20-40 /im (the pellet diameter is about two nozzle 
diameters) has been observed. Deviations from the mean 
pellet diameter are below 1% (10%) over periods of few 
seconds (hours). The average velocity of 30 /im pellets 
amounts to ~70 m/s. The radial displacement of the pel- 
lets from their nominal flight path in the dummy chamber 
of about ±200 /im has been extrapolated from the angu- 
lar pellet distributions measured behind the 1 st sluice. 
This value is dominated by our experimental resolution. 
The small pellet-beam divergence is a consequence of the 
extremely regular jet breakup in the TPC (Fig. [3]). 

If used as an internal target at a storage ring, our 
device minimizes unwanted gas loads to the accelera- 
tor and allows for an effective target thickness of a few 
10 15 atoms/cm 2 . Since the tube connecting the genera- 
tor with the interaction zone is rather narrow, detectors 
can be placed close to the interaction point in a nearly 4tt 
configuration. Such a geometry is also advantageous for 
placing collector mirrors in laser-plasma UV sources Q. 
In addition, the large separation of the hot and radiat- 
ing plasma from the jet generator makes such a set-up 
especially favorable. 

In this letter we present first data on the mono-disperse 
breakup of cryogenic jets. We find that the jets are sig- 
nificantly more resistant to breakup than predicted by 
Rayleigh's theory. We suggest that a reason for this dis- 
crepancy is the influence of evaporation effects. This 
interpretation is supported by the observation of non- 
axially symmetric jets. This new jet mode — for which 
so far not even a rudimentary theoretical explanation has 
been formulated — sets in when one decreases the ambi- 
ent pressure very close to the triple-point values. For a 
complete theoretical delineation of drop formation pro- 
cesses from thin cryogenic jets, further systematic data, 
like on the X dependence of 7/70 and the pressure de- 
pendence of Ljet, are required. 
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